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Q ' Abstract 

^ . In this paper a neutron star with an inner core which undergoes 

bJQl a phase transition, which is characterized by conformal degrees of 

freedom on the phase boundary, is considered. Typical cases of such 
CN ' a phase transition are e.g. quantum Hall effect, superconductivity 

, and superfluidity. Assuming the mechanical stability of this system 

Q>^ I the effects induced by the conformal degrees of freedom on the phase 

boundary will be analyzed. We will see that the inclusion of conformal 
degrees of freedom is not always consistent with the staticity of the 



I phase boundary. Indeed also in the case of mechanical equilibrium 

■ there may be the tendency of one phase to swallow the other. Such a 

. shift of the phase boundary would not imply any compression or de- 

compression of the core. By solving the Israel junction conditions for 
the conformal matter, we have found the range of physical parameters 
which can guarantee a stable equilibrium of the phase boundary of 
5^ I the neutron star. The relevant parameters turn out to be not only the 

density difference but also the difference of the slope of the density 
profiles of the two phases. The values of the parameters which guar- 
antee the stability turn out to be in a phenomenologically reasonable 
range. For the parameter values where the the phase boundary tends 
to move, a possible astrophysical consequence related to sudden small 
changes of the moment of inertia of the star is briefly discussed. 



Keywords: Junction conditions, Neutron Stars, Conformal Boundary De- 
grees of Freedom. 
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1 Introduction 

One of the most intriguing features of neutron stars is the possibihty of having 
an inner core which undergoes a (colour) superconductivity and/or superflu- 
idity phase transition due to the extremely high pressure and density. This 
fact was first pointed out by Migdal [1] (there is a huge amount of literature 
on this subject with a little hope of providing one with a complete list of 
references: for two detailed reviews, see |2] [5] and references therein; for 
general relativistic formalisms suitable for dealing with neutron stars within 
these scenarios see [1] [5] and references therein). Such phase transitions in- 
side neutron stars could lead to interesting observational effects relevant in 
cosmology [6] and through the quasi normal modes of the stars as well as the 
related gravitational radiation (see, for instance, ^ detailed reviews on 
the subjects are |9] [10] [H]). Recently, there has been also pointed out the 
intriguing possibility to have a quantum hall phase of gluons and quarks in 
the inner core of a neutron star (see, in particular, [12] [13]). The discovery 
of magnetars [T?], highly magnetized neutron stars whose magnetic fields 
can reach lO^^G in the core, makes the arising of Quantum Hall features in 
the inner core of a neutron star not unlikely. All the here described types 
of phase transitions have in common to be characterized by the presence of 
conformal degrees of freedom on the phase boundary predicted by QFT0 see, 
for instance, p3J [fS])- 

The main goal of this paper is to study under which conditions the phase 
boundary inside a neutron star, which undergoes such a phase transition with 
conformal boundary degrees of freedom, can be static in general relativity, 

^The fact that the boundary degrees of freedom are conformal can be seen "heuris- 
ticahy" in the case of a superconducting phase. Due to the Meissner effect, the current 
must circulate on the phase boundary of the superconductive phase. Phenomenologically, 
there is no dissipation, this means that there is no physical scale over which the boundary 
excitations can dissipate energy. The lack of a characteristic scale for the boundary theory 
is related to conformal symmetry. In the case of Quantum Hall Effect, this can be easily 
seen from the fact that it is well described by a Chern-Simons action in the bulk. Thus, 
the theory induced on the boundary is the Wess-Zumino-Witten action which also has 
conformal symmetry. 
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assuming mechanical equilibrium of the systenH. Indeed also if the system 
is mechanically stable i.e. there is no pressure inducing the collapse or ex- 
pansion of the core, there can be another sort of instability: there may be 
the tendency of one phase to prevail over the other near the boundary. This 
would imply that the boundary gets shifted. This situation is analogous to a 
bubble chamber where a small perturbation induces a phase transition. To 
study this phenomenon one must take into account the contribution of the 
nontrivial traceless boundary stress tensor to the Einstein field equations. 
This problem is equivalent to solving Israel's junction conditions. It is nec- 
essary to assume a dynamical phase boundary and to check under which 
conditions such a configuration is a solution of the Einstein field equations 
(by solving the junction conditions). It turns out that the equation of mo- 
tion for the dynamical phase boundary is equivalent to the dynamics of a 
classical point particle in an effective potential. There exist a static stable 
equilibrium configuration if the effective potential has a local minimum for 
a suitable negative value of the potential. It is shown that the existence 
of a local minimum is regulated both by the difference in the mass density 
between the two phases and by difference of the slope of the two density 
profiles. A third parameter which determines the form of the potential is the 
energy of the shell. The non trivial effects of conformal boundary degrees of 
freedom have been pointed out, in the context of black hole physics, in [19] 
in which they are related to the arising (after the event horizon is formed) 
of Quantum Hall features [20] [U] related to the strong attractive nature of 
the gravitational field acting on Fermions inside a collapsed neutron star. 

The structure of the paper is the following: in Section 2, the assumptions 
of the present paper are explained. In Section 3 the junction conditions 
derived from the Einstein equations for a neutron star with a phase boundary 
are solved. In section 4 the range of the physical parameters, characterizing 
our configuration, for which the phase boundary is in stable equilibrium is 
determined. In section 5 an astrophysical implication is discussed. Section 6 
the conclusions and perspectives are drawn. 

^ The analysis of the mechanical stability of neutron stars has been performed in a 
huge number of papers. To provide one with a complete list of references is a completely 
hopeless task; representative papers are, for instance, [T^ [IS] |3j and references therein) 
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2 The standard approximations in a neutron 
star 



Inside a neutron star, at the energy scale of the standard model of particle 
physics (up to TeV) the collapsing neutron star can be described very well by 
QFT and classical general relativity. The quantum dynamics of the neutrons 
and of the quarks living inside the neutron star is much faster than the 
dynamics of the gravitational field. This implies that one can compute the 
equation(s) of state of the Fermions as usual and then use such equation(s) to 
solve the Einstein equations in which the source is described by the equation 
of state itself. The success of this theory initiated by Landau, Chandrasekhar, 
Tolman, Oppenheimer, Volkoff, Snyder (and many others) tells that such an 
adiabatic approximation is excellent. The typical mass M^s, radius Rns, 
density p^vs and the baryon number N]srs a neutron star are 

M^s ~ 10'=^^?, Rns ~ (lO' ^ 10^) cm, 
PNs ~ {10'' ^ 10'') g/cm^, N^s ^ 10'\ 

A neutron star pulsation can be thought of as a linearized solution of 
the Einstein equations around a background describing a static star (see, 
for the spherically symmetric case, Eq. ([1])) with a "quasi periodic" time 
dependence: that is, the perturbation h^^, of the background metric Eq. ([T]) 
has the following form 

V t) = exp (^i (^u + ?j U'^J (^) 

where the imaginary part of the frequency ^ is due to the damping by emis- 
sion of gravitational and electromagnetic waves (see, for instance, [9]). The 
standard order of magnitude of the frequency of a neutron star pulsation is 

19, 85Hz < < 12, 84kHz 

while the damping has a wider range of possible values varying from low 
damping (in which r can be of the order of years) to high damping (in which 
T is of the order of millisecond). 
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3 Junction conditions and equation of motion 
of the phase boundary 

In this section the junction conditions for the considered neutron star are 
found and the physical parameters that determine the existence of stable 
equilibrium configurations are found. We will set 

h = 1, c = 1, 

while keeping explicitly the Newton constant. 

In the present case, one would like to describe a star in which the inner 
core underwent a gapped phase transition and therefore conformal boundary 
degrees of freedom localized at the separation between the two phases (pre- 
dicted by QFT [U] [IS]) have also to be included. We will therefore assume 
that there are two phases, each in mechanical equilibrium, namely the inte- 
rior region which represents the exotic (superfluid, superconducting quantum 
Hall) phase and the exterior region which is the normal phase. Each phase 
is approximated by a perfect fluid in a static configuration. The two regions 
are characterised by very different equations of state. 

We shall see that, even assuming static equilibrium (without any compres- 
sion or decompression of the core), the conformally invariant matter living 
on the phase boundary may not be consistent with a static junction. In 
this case, since we assume that the conformal invariance is a defining char- 
acteristic of the phase transition, there must be some tendency of the phase 
boundary to move. Such a change in the position of the phase boundary is 
not by collapse of matter but by a kind of "creeping" phase transition, where 
one phase tends to swallow the other. 

The metric describing a spherically symmetric neutron star can be parametrized 
as follows 

ds^ = - exp{2u)dt^ + exp{2X)dr^ + r^dQ^ , (1) 
u = z/(r), A = A(r) 

where dil"^ is the line element of a unit spher^. The functions u and A are 
determined by the Einstein equations. If the source is a perfect fluid one gets 

■^The description of a rotating neutron star is much more difficult. However, at least 
for slowly rotating neutron stars, one can argue that the analytic results derived here in 
the spherically symmetric case do not change qualitatively as long as the "rotation" can 
be dealt as a perturbation [10] [26]. 
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the standard Tolman-Oppenheimer-Volkoff-Snyder equations: 



exp(2A) = (l - 2G^^^ , m(r) = 4vr ^ x^p{x)dx (2) 
d,P = -(p + p)dru, g,z/ = exp(2A) (3) 

where p and P are the density and the pressure of the perfect fluid. 

One has to find under which conditions such a configuration solves the 
Einstein field equations. Given the static perfect fluid solutions in the interior 
and exterior, this problem is equivalent to solving Israel's junction conditions 
|21j . This technique allows the construction of an exact solution of the Ein- 
stein equations by matching two different solutions provided that the metric 
is continuous and the discontinuity of the extrinsic curvature is compensated 
by a suitable energy momentum tensor S'^ describing the boundary degrees 
of freedom of the inner phase: in the present case in which the inner phase 
transition is gapped one should only allow traceless S^. Let S be a timelike 
hypersurface of codimension one, on which the matching has to be performed. 
Let n'^ be the (spacelike) unit normal to E, and let h^'^ = g^^ — n^n" be the 
metric induced on S. The matching conditions at S are: that the induced 
metric be continuous; that the jump in the extrinsic curvature is related to 
the intrinsic stress tensor localized on S by [21] 



7^; = -87rG(5^:-^5>5), (4) 



where 7^,^ is defined as 

7^,, = lim [K^^{r] = +e) - K^,^{r] = -e)] 



+0 



7] being the arc length measured along the geodesic orthogonal to S. The 
extrinsic curvature Kf^^, of S is given by 

Ko,p = h>'^h';V,n,. (5) 

In the present case, the matching is at a timelike hypersurface r = r*(r) 
(let us call such a surface ^{t)) which describes the time evolution of a thin 
spherical shell where the phase boundary occurs. Namely, at r = r* a gapped 
transition occurs (such as superconductivity, or Quantum Hall state). 
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S(r) 



3.1 The phase boundary 

Because of the spherical symmetry, it can be assumed that the spatial sections 
of the timelike matching hypersurface E(^) are isomorphic to the two sphere 
so that there exist a coordinates system in which the induced metric ds^ ' 
and 5*^ respectively read: 

Sf"" = a{T) (UW) - C{t) {h'"' + WU") , (6) 

where h'^" is the metric induced on S(^), n'^ is the (spacclike) normal to S(t-), 
T is the arc length measured along the timelike geodesic belonging to S(t-), 
U^^ is the normalized four velocity of E(^), a is the surface energy density of 
^(r)) C (''") is the surface tension and the conservation of implies 

r 

with r(T) being the proper circumferential radius of the domain wall E(^).The 
requirement of a traceless S^" (otherwise it would not correspond to the 
classical description of boundary gapless degrees of freedom) leads to 

a + 2C = (7) 

which combined with the above conservation equation implies 

a=^, ao>0 (8) 

where (Tq is an integration constant, with the dimension of an energy times a 
length, which depends on the microscopic model and which will be approxi- 
mated by a characteristic energy scale of the gapped phase transition times 
the length scale of the domain wall. 

To get an intuition on the meaning of ctq let us consider the case in which 
in the inner core of the neutron star there is a Quantum Hall phase. In such 
a case, the only physical meaningful energy scale of the model is the energy 
gap Egap, thus one can expect that 



ctq ~ I Egap d{i 
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where r* is a typical value of the radius of the inner core, the index i labels 
the Landau level^, d{i) is the degeneracy of the i-th. level, and nmax is the last 
fully filled Landau level (so that (Tq/ (r*Egap) represents the number of parti- 
cles living in the Hall phase). In the case of a supeconductivity/superfluidity 
phase inside the inner core also one should expect a close relation between 
(To/r* and the energy gap times a number of the order of magnitude of the 
number of particles living in the gapped phase. 

The "equation of motion of the domain wall" , that is the equation which 
determines the evolution of t'(t), can be deduced from the matching condition 
dlj). In particular, the important equation is the angular components SI of 
Eq. (HD 

7i = -SnGSi (9) 
One can see that the rr component of the matching equations 

7; = -8nGS; (10) 



is not independent of the angular ones provided the consistency condition [2T| 
^^\u = ^ n'^T^j,y(ext) — n'^T^j,(int) = for the conservation of energy on 
S is taken into account: one can easily see that the rr component of the 
matching equations only constraints the discontinuity of the first derivative 
of the metric function u. Note that the pressure and density profiles are 
assumed to be static. 

The standard procedure to compute Yi is as follows: the unit normal to 
the junction surface is 



= ± (^-exp(i^ + A)5,r,exp(2A)-\/(a,r)' + exp(-2A),0,0j (11) 

where the ± ambiguity depends on whether one considers the inner or the 
outer directed unit normal. A straightforward calculation using ([5]) shows 
that the angular components of the extrinsic curvature is 



^In contrast with the ordinary Quantum HaU devices in which one often can assume 
that only one Landau Level is fully filled, in the present case (because of the huge density 
inside the inner core of a neutron star) one should expect that many Landau level may be 
fully filled. 
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while the rr-component of the extrinsic curvature is ( [22] ) 

v' {{drrf + exp(-2A)) + (9» + A' {drvf 



Eventually, Eq. ([9]) reads 



{d^rf + exp(-2A) 



(13) 



(i^(.nt)-i^(ext)) = ^.^ 



(14) 



where there is no sum on the repeated indices in this formula, K\{iv£) and 
Kl{ext) stand for the extrinsic curvature computed on the inner and on the 
outer side of the junction. This equation can be written as an ordinary first 
order equation for r(r) bringing on the right hand side Kl{ext) and then 
squaring: 



(drr) 



(exp(-2Ai„t) - exp(-2Aext)) 



K 



1 2 



2K ■ 2r2 

which, using ([2]) can be written: 

{drvf + V{r) = , 

2Gmextir 



V{r) 



-1 + 



(15) 

(16) 



where Xint and Xext stand for the metric functions of the inner and of the 
outer side of the junction and niextij) is the mass function ([2]) of the exterior 
part. It is worth to note an interesting feature of the present description: the 
dynamics of r(r) given by equation (IT^ only depends on the metric function 
exp(— 2A) which in turns only depends on the density (see Eq. ([2])). 

The pressure inside and outside of the shell is given by the Tolman-Volkov- 
Oppenheimer-Snyder equations 



drP 



(p + P) 1-2G 



(m + Anr^P) 



(17) 



where m and p on each side are functions of r. Therefore we get an equation 
of the form drPextO-ir) + h{r) + c{r) Pextij") + d{r)P^^^{r) = (and similarly 
for Pint), where a, b, c and d are polynomial in r, has a solution (which may 
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be found numerically) with one constant of integration. For a consistent 
solution, we must take into account equation (ITU]) . To see the existence of a 
solution it is sufficient to consider the static case at the equilibrium radius 



Teq, in which the equation reduces to: 



eq 



2m{req) 



-I ext 



int 



^eq 



(18) 



where 



1 ext 
lint 



means the jump in this quantity across the shell. Generally, 



we can see this as an equation constraining the two integration constants for 
Pint and Pext in terms of K. It is enough to know that there is a value for 
the pressure, but the precise value shall not be needed to solve the stability 
problem, since the potential depends explicitly on p and not P. 

One sees also that the equation of motion f|T5l) for the shell r(r) is formally 
analogous to the equation of motion of a classical point particle where V{r) 
takes the role of an effective potential. A necessary condition to have a static 
phase boundary is that V{r) must have a local minimum. The exact form of 
the potential depends now on the functional form of the density profiles of 
the two phases. 

It will be now explicitly shown that in a suitable range of parameters such 
a local minimum exist. 



4 Density profiles and (non)staticity of the 
phase boundary 

To find the possible relative minima (if any) of the potential V{r) one has 
to write explicitly Xint and Xext as functions of r and of phenomenological 
parameters characterizing the density of the star. At first glance, to do 
this one should solve the Tolman-Oppenheimer-Volkoff equations inside and 
outside the inner core. This is a very hard task since to do this one should 
approximate the equation of state of strongly interacting quarks, gluons, 
neutrons and so on at very high densities and pressures. In fact, this can 

^An analogous situation possibility has been studied [T^ in the context of black hole 
physics, where the conformal degrees of freedom led to the existence of a local minimum 
for the effective potential (In [23] , where a boundary stress tensor of a cosmological term 
is considered, there was no local minimum). 
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be avoided: the reason is that the effective potential only depends on the 
metric functions Aj„i and Xext (that is, on grr)- Such function is expressed 
in terms of the density profile: if one is able to determine the dependence of 
the density on the radial coordinate r the effective potential can be written 
explicitly. Remarkably enough, in many sound models of neutron stars it is 
possible to determine the density profile explicitly even if the knowledge of 
the equation of state (that is, the functional relation between density and 
pressure) is not perfect (see, for instance, [5] [10] [23] and references therein). 
On the basis of the above references, it will be here assumed that the density 
of the neutron star has the following form: 

p.Ur) = ^-^r", r<r(r) (19) 
PeAr) = r>r(r), (20) 



where 



n = l, Ainu Ae.t ~ (10" ~ 10^^) g/cm\ B^nu Be.t ~ (lO" ~ 10^^) g/cm\ 

Of course, even if it is not explicitly present, the pressure is needed to deter- 
mine the values of the above density profiles. The choice of a linear decreasing 
is general enough since in the present case the important region is the one 
near the boundary phase so that a linear decreasing captures the main qual- 
itative features (for instance, if one would take a quadratic decreasing (that 
is, n = 2 in Eqs. f|T9l) and fl20l) ) as in [21] the main conclusions of the present 
paper would not change). On the other hand, interesting phenomenological 
conclusions can be reached in the cases in which further terms are added in 
the density profile. Such consequences will be discussed in a moment. There- 
fore, one can assume that both in the inner core (r < r(r)) and outside the 
core ( r > r(r)) the density decreases linearly. With the above choice of the 
parameters the effective potential (fT6l) reads 

V{r) = -(2G (^r^ - - 1 + 



1 

+ 



^4: 



Q I Ag^xt ~ Aint \ ^ f Bgxt ~ Bint \ ^7 I K 



3K \ AK 2 



(21) 



To proceed, let us adopt the geometric unit of measure in which any quantity 
is expressed in powers of cm (and the Newton constant is 1): 

1 MeV ^ 10"^^ cm, G = l. 
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Typical values are: 
having assumed 

Egap ^ MeV, ^ 102^ r* ^ cm 

where N is the order of magnitude of the number of particles living in the 
inner phase and r* is the typical radius of the inner core. 



4.1 Stable range 

From a numerical analysis, it can be seen that in order for there to be a 
minimum of the potential we must have the following relations 

AAAB > 0, 

a further condition is that 

{drrf + V{r) =0 V{r) < 0. (22) 

When the above condition is not fulfilled, the Israel equations cannot be 
satisfied at all (the matching would be impossible). The numerical graphs 
also suggest that for a given AB there is a minimum value of AA below 
which solutions do not exist. For a core of radius ~ 1 cm this minimum is 
AAmin ~ AB/2 as shown in graph [2l More generally, one can check the 
following: 

AM = — - > Kr' . (23) 

3 4 - ^ ^ 

This is natural (by analogy with special relativity) since Kr^ is the rest mass 

of the shell and Am is a measure of the relativistic mass. If this condition is 

not fulfilled, there is either no minimum of the potential, or the minimum is 

positive, meaning no stable solution (see figure 2 ). 

In the range of parameters we focus on (see graph [1]) a stable core exists: 

the equilibrium radius r^g and the frequency Upui are 

(24) 



eq 


^0.53 cm. 


dlv\ 

' ^min 


= 146cm ^ , 


eg 


^1.32 cm. 


dlv\ 

' ' mm 


= 23.1cm"2 


eq 


^3.32 cm. 


dlv\ 

1 T - 

' ' mm 


= 3.67cm~^ 


eq 


^5.26 cm. 


dlv\ 

' ' mm 


= 1.46cm~^ 


eq 


?«8.34 cm. 


dlv\ 

' ^min 


= 0.582cm- 
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Figure 1: In this figure there is a graph of the effective potential for varying 
10"^cm"^ < AB/iK < lOOcm"^ and suitable values of AA/3K such that the 
potential has a minimum for small negative V. Stable solutions with radius of 
order 1 — 10 cm are shown to exist. 



So the frequencies are of order of GHz. 

The above frequency of pulsation would correspond to a periodic motion 
of r(r) between a maximum and a minimum radius oscillating around the 
above equilibrium radiu^. The equilibrium radius considered here is of the 
order of centimeters (which is in good agreement with the QFT estimates in 
the case of the Quantum Hall Effects [13] [12]; such a value for the radius of 
the inner core could also be compatible with the presence of a superconduc- 
tivity/superfuidity phase [2] [3]). 

4.2 Nonstatic range 

The numerical analysis shows that if AA and AB do not have the same sign 
the configuration is unstable: that is, the effective potential (!2T!) has only 
one local maxima and the phase boundary will tend to move. It is important 
to stress that this result is a purely general relativistic effect related to the 
conformal boundary degrees of freedom. Such an effect would be overlooked 

^Tlie actual possibility that such an exact pulsation mode could be realized in practice 
in a typical neutron star will be discussed in the next section. 
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Figure 2: In this figure there is a graph of the effective potential for fixed 
l^B/AK = and varying values of /S.A/'iK (see plot legend above). Stable 

solutions exist for IS.A/?>K greater than approximately 65cm~^. 



if one neglects the possible dynamics of the phase boundary. If one construct 
a neutron star by matching the inner core and the outer part at a fixed 
radius r* = const (neglecting its possible dependence on r) and without the 
inclusion of the boundary degrees of freedom then it is not possible to disclose 
the effective potential. 



4.3 Multiple equilibria? 

As it has been already discussed, a density which decreases linearly with the 
radius is reasonable and general. On the other hand, it may be well possible 
(see [5] [To] and references therein) that the density profile is corrected by 
further terms as for instance 

P^nt{r) = ^-^r + Q^tr'', r < r(r) (25) 





Bint 


Air 


Att 


Aext 


Bext 


An 


Att 


1. 





Pextir) = ^-^r + D,,,r-, r > r(r), (26) 

(27) 
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At a first glance, one could expect that near the phase boundary only the 
first two terms in Eqs. (1251) and fl26p are important. However, if due to the 
extreme conditions inside a neutron star, such further terms have to be taken 
into account, then interesting phenomenological scenarios open up. In these 
cases, the effective potential may have multiple local minima and the valleys 
would be generically asymmetric: namely the valley corresponding to the 
larger radius of equilibrium could be more deep than the valley corresponding 
to the smaller radius of equilibrium. In principle, a phase transition could 
occur in which the phase boundary jumps from one valley to the other. 

5 Astrophysical Implication 

From the explicit form of the effective potential (12T|) it is clear that, in the 
range of parameters in which the effective potential has a local minimum, it 
is consistent to assume a static phase boundary. In this range of parameters 
also small oscillations of the phase boundary could be possible. However it is 
difficult to say what would be the sources of dissipation which would damp 
the oscillations. 

It is a well known fact since the Migdal paper [1] that the presence of a 
superfiuid phase inside a neutron star affects its moment of intertia. In the 
range of parameters in which the phase boundary tends to move one phase 
tends to swallow the other. In the case in which the inner phase expands, 
there is a corresponding decrease of the moment of intertia. Assuming that 
our results also hold for the rotating case (which is a reasonable assumption 
at least in the slow rotating case [lO] [26]) this would lead to interesting 
observable effects: namely a sudden small increase of the angular velocity 
because of the angular momentum conservation. Analogously, in the case in 
which the inner phase shrinks, there is a corresponding decrease of the angu- 
lar velocity. Indeed both types of phenomena are of interest in astrophysics. 

6 Conclusions and perspectives 

We have studied a model of a neutron star with an inner core which under- 
goes a phase transition. It has been shown that, even assuming mechanical 
equilibrium, it is not always consistent to assume a static phase boundary 
once conformal boundary degrees of freedom are taken into account. 
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In order to study the staticity of such a phase boundary one must take 
into account the general relativistic effects of these boundary degrees of free- 
dom by including a nontrivial stress tensor on the junction. To the best of 
the authors knowledge, such a consistency analysis with conformal boundary 
degrees of freedom has not been performed previously. The astrophysical 
consequences of the non-static regime, related to sudden changes of the mo- 
ment of inertia of the star, are worth further investigating. 
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